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Abstract 

In this work we investigate the energy gap between the ground state and the first 
excited state in a model of two single-mode Bose-Einstein condensates coupled via 
Josephson tunneling. The energy gap is never zero when the tunneling interaction 
is non-zero. The gap exhibits no local minimum below a threshold coupling which 
separates a delocalised phase from a self-trapping phase that occurs in the absence 
of the external potential. Above this threshold point one minimum occurs close 
to the Josephson regime, and a set of minima and maxima appear in the Fock 
regime. Expressions for the position of these minima and maxima are obtained. The 
connection between these minima and maxima and the dynamics for the expectation 
value of the relative number of particles is analysed in detail. We find that the 
dynamics of the system changes as the coupling crosses these points. 
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1 Introduction 



Research into Bose-Einstein condensates (BECs) continues to attract considerable atten- 
tion from both theoretical and experimental sides due to its contribution to the compre- 
hension of physical phenomena emerging from mesoscopic systems. Recently, with the 
modernisation of experimental techniques, BECs have been produced and investigated 
by several research groups using magnetic traps with different kinds of confining geome- 
tries [1,2] and chip devices [3,4]. From the theoretical point of view, the model of two 
single-mode BECs coupled via Josephson tunneling (see eq. [T]in the next section and a 
comprehensive review in [5]) has been a successful model in further understanding the 
experimental realization of BECs in dilute alkali gases [6]. The success of this model is 
based on two main reasons: first is that despite its apparent simplicity, the model captures 
the physics behind the macroscopic quantum phenomena that emerge from these systems; 
second is that it is also applicable in other studies in mesoscopic solid state Josephson 
junctions [5,7,8], non-linear optics [9], statistical physics of spin systems [10,11] and in 
nuclear physics [12,13]. A theoretically important feature of this model is it can be formu- 
lated through the Quantum Inverse Scattering Method establishing solvability by Bethe 
Ansatz methods [14,15]. 

Many aspects of this model have already been investigated, such as the dynamics of 
the population imbalance [16,17,19], the imbalance fluctuation [18,19] and entanglement 
[9, 19-22] amongst others [23-25]. However an understanding of the properties of the 
energy gap as a function of the external potential, and an investigation of the dynamical 
consequences, appears lacking. This study is of interest since it could feasibly be tested 
in an experimental setting, particularly in cases where the external potential is the only 
tunable interaction, as may be expected in solid state systems. Here, we will expand on 
this aspect of the model by undertaking a detailed analysis of the energy gap behaviour 
across all coupling regimes. Specifically we will show that the energy gap function presents 
no local minimum below a threshold coupling, defined below, while above this threshold 
point one minimum occurs close to the Josephson regime and a set of minima and maxima 
appear in the Fock regime. The connection between these minima and maxima and the 
dynamics for the relative number of particles is investigated. We find that these minima 
and maxima act as a transition point separating entirely different dynamical behaviours. 

The paper is organised in the following way. In section II we present the model and 
discuss the energy gap for different choices of the tunneling coupling and different regimes. 
In section HI we discuss the correspondence between the minima and maxima that occur 
in the energy gap function and the temporal evolution for the relative number of particles 
in different regimes. In section IV we compare the quantum results with those obtained 
from a classical analysis. Section V is reserved for our conclusions. 

2 Energy gap 

The model we will study, based on the two-mode approximation as described in [5] (see 
also references therein), is described by the following Hamiltonian: 
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where al,ai {i = 1,2) denote the single-particle creation and annihilation operators as- 
sociated with two bosonic modes and A^^i = a\ai and N2 = a2'^2 are the corresponding 
number operators. The quantity = A^^i -|- A^2 is the total boson number operator and is 
conserved. The coupling k provides the strength of the interaction between the bosons, 
A/i is the external potential and Sj is the coupling for the tunneling. We mention that 
throughout we will take k > 0. The change Sj —Sj corresponds to the unitary trans- 
formation CLl — > Oil, CL2 — ^ — (l2i 

while A/i — ^ — AyU corresponds to ai ^ a2- Therefore we 
will restrict our analysis to the case of Sj, Afi > 0. It is useful to divide the parameter 
space into a number of regimes; viz. Rabi {k/£j < 4:/N), Josephson < k/£j << N) 
and Fock {k/Sj » N). In the Rabi and Josephson regimes, coherent superposition of 
the two condensates is expected whereas in the Fock regime the two condensates may 
be essentially localised. Note that we adopt the convention ( cf. [26]) to specifically de- 
fine the boundary between the Rabi and Josephson regimes by the threshold coupling 
k/Sj = 4/N. This threshold was found in [16] as the transition coupling between macro- 
scopic self-trapping {k/Sj > 4/N) and delocalisation k/£j < 4/N) for the case A/i = 0. 
(see also [19,27,28]). The presence of a non-zero external potential changes the coupling 
ratio at which this transition occurs, as we will show. However, it is useful to still refer to 
k/£.j = 4/N as the threshold coupling, irrespective of the value of the external potential, 
for other reasons which we will discuss. 

Our main aim here is to investigate the behaviour of the energy gap for the Hamiltonian 
(^. In the Fock regime, we can preliminarily study this problem by putting the tunneling 
interaction term of the Hamiltonian to zero, i.e. 8j = 0. In this case all energy levels 
can be computed exactly. For the case of an even number of particles this yields 

^ 9 . N N N N 

Ep = -v" + pA/i, p = --, -- + 1, -1, 0, 1, 1, -. (2) 

In the absence of the external potential (A/t = 0) the ground state is unique {Eq = 0) and 
all other energies are double-degenerate, since E_p = Bp. By switching on the external 
potential A/i, these degeneracies are broken. 

Fixing the total boson number N and the interaction strength k and allowing the 
external potential A/i to vary, expression (j2I) indicates two important facts: i) there are 
no degenerate energy levels, only level crossings; ii) all energies are functions of the value 
of the external potential value, except for Eq = 0, and all energies of the type E^p and 
Eq correspond to the ground state in some region, defined by the strength of the external 
potential A/t. These regions are determined by the level crossings. Results are depicted 
in Fig. 1 for the case iV = 6, A; = 6. 

The energy gap, in each region, is given by 

Ay+i,,-) = - E.j , A/i G [kj, k{j + 1/2)] 

A(,_i,,) = , A/iG [A;(j-l/2),A;j] (3) 

A(Ar/2-l,Ar/2) = E^(^n/2~1) - E_n/2 , A^>k{N-l)/2. 

Here A(jj) (A^-^j)) denotes the energy gap before (after) the level crossing between the 
energy levels E^i and E^j. Notice that at each different region there is either a different 
ground state or a different first excitation. A similar structure has also been found in [29] 
by fixing all coupling parameters and varying the total number of particles. 
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Figure 1: Energy levels versus the external potential for = 6 and k = 6. The level 
crossings involving the lowest energy level represent changes in the ground state structure. 
The energy gap is also determined by level crossings involving the first excited state. 



Considering the case of an odd total number of particles A^, eq. 1^ is still correct, 
however we have to eliminate the energy level Eq and change the allowed values of the 
index p to p = -N/2, -N/2 - 1, -1/2, 1/2, N/2 - 1, N/2. Here, at A/i = 0, all 
energies are two-fold degenerate. By switching on the external potential A/i, all degen- 
eracies are again broken and level crossings appear. Comparing the two cases, one can 
see that the gap behaviour of the model with external potential A/i and an odd number 
of particles N is the same as the profile for a system with + 1 particles and external 
potential A/i + k/2. 

Now we consider the situation where the coupling Sj between the two condensates is 
switched on and we study its effect using the numerical diagonalisation of the Hamiltonian 
(^. In terms of the energy gap, the existence of the threshold coupling is suggested by 
the curve of the energy gap as a function of the external potential. Above the threshold 
[k/Sj > 4:/N), one or more local minima appear while no local minimum occurs when 
k/£j < i/N. Let us first focus our attention on the analysis of the energy gap in the 
region above the threshold point, in particular, in the region where k/Sj > 1. In Fig. |21 
we plot the energy gap A versus the external potential A/x with different values of the 
coupling parameter Sj for the cases A^ = 6, 7. 



We can see clearly the effect produced by the tunneling term Sj. Before switching on 
this coupling, in the extreme Fock regime {Sj = 0), the energy gap curve is not differen- 
tiable at each minima and maxima. Here the maxima and minima are determined by level 
crossing and, consequently, either a change in the ground state or the first excited state. 
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Figure 2: Energy gap A versus the external potential A/x for different choices of the 
coupling parameter £j = 0, 0.2, 1, 5, 6. On the left, N = 6, and on the right, N = 7. In 
the extreme Fock regime {Sj = ) the minima and maxima represent the level crossings as 
depicted in Fig. 1. For small values of Sj the difference between two consecutive minima 
or maxima is constant and equal to k. As Sj increases just one minimum survives. 



In this particular situation we can exactly calculate the energy gap value in all regions 
using equations (j21 El)- It is interesting to observe that at each minimum (respectively 
maximum) the same value for the energy gap is found. 

Turning on the coupling term Sj the energy gap function becomes smooth. An im- 
portant fact here is that weak tunneling does not considerably change the position of the 
minima and maxima in relation to the Sj = case. It does change, however, the value 
of the energy gap at each minimum and maximum. This effect is strongest in the first 
region (A/x G [0, fc/2]), reducing gradually at each subsequent region, with the last min- 
imum corresponding to the lowest value of the gap. Increasing the value of the coupling 
term Sj to take the model out of the Fock regime leads to a suppression in the number 
of minima and maxima until just one minimum survives. 

We now look closer at the set of minima and maxima that appear in the Fock regime. 
The number of these minima and maxima depends on the total number of particles A^. 
If A^ is even, the number of maxima and minima is the same, A^/2, while if A^ is odd the 
number of minima is (A^ + l)/2 and the number of maxima is (A^ — l)/2. However, we 
find that the difference between two consecutive minima (or maxima) is always constant 
and equal to k independent of the total number of particles. This allows us to examine 
a general expression for the position of each minimum and maximum of the energy gap 
that occurs with respect to the external potential, which we find fits very well with 

Afi, = ^^k + ^, 1 = 1,2,3,. ..,N. (4) 
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Above, odd choices for the index / (/ = 1, 3, 5, 7, ...) generate the minima, while even 
choices (/ = 2,4,6,8,...) the maxima. The particular choice / = 1 corresponds to the 
lowest minimum, while / = 3 the second lowest minimum, and so on. With the exception 
of the lowest minimum, all other minima and maxima disappear as the ratio £j increases. 
This is illustrated in Fig. |21 for = 5, 6 corresponding to the situation where we approach 
the Josephson regime after leaving the Fock regime. For k/£j > 1 (Josephson and Fock 
regimes), we find that the position of the minimum can still be determined by the previous 
expression Q using 1 = 1. 

In the interval k/Sj G (4/A^, 1] the energy gap function still has a local minimum, 
however the position at which this occurs is not longer given by and we have not been 
able to find a simple expression for it. The value of Afi at which the minimum occurs 
approaches zero as the ratio k/Sj approaches the threshold coupling. Further decreasing 
the ratio k/Sj past the threshold point toward the Rabi regime, we find a new scenario: 
the energy gap is now a monotonic function for all values of A/i. This behaviour is 
illustrated in Fig. |2l where we plot the energy gap function versus the external potential 
A/i in the region k/Sj < 1 (Josephson and Rabi regimes, including the threshold point) 
for N = 40 particles. In these cases differences between even and odd particle numbers 
are not significant as the strength of the tunneling term dominates any parity effects due 
the external potential. 



N = 40 




AM./k A(j./k 

Figure 3: Energy gap versus the external potential for = 40 and different choices of 
the parameter ratio £j/k = 1, 2, 5, 9, 10, 11, 15, 20. On the left (right) side we display the 
points below (above) the threshold point Sj/k = 10. On the left hand side we can clearly 
see the existence of a minimal energy gap, while on the right hand side the energy gap is 
a monotonic function of the external potential. 
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We emphasise that (jlj) was obtained by numerical fitting. It is useful to see how this 
result compares with the predictions of perturbation theory. It was argued in [14] that 
for the weak tunneling regime perturbation results are only valid for the Fock regime 
Sj << k/N. A standard calculation using degenerate perturbation theory shows that 
to first order there is no correction to the values of the applied potential which give the 
minimal energy gaps; i.e. 

- / 

^A*pert = — 2 — ^' ^ = 1, 3, 5, 7, . . . 

To make comparison between (0]) and this perturbative result we note that for the Fock 
regime 

N-r 



-k < Afic 



N -I k 
N-l 1 \ 

N-l ^ [N-l 1 \ 
^A:<A/..<^^ + -j k. 

Thus the formula Q agrees with the perturbative result up to a small correction of 
order A^~^. It is important to stress however our numerical result ^ applies also in the 
Josephson regime where the perturbative result is not valid. 

It was also argued in [14] that perturbation results for strong tunneling are only valid 
for £j» kN . Our numerical results show that the existence of minima for the gap 
only occurs for couplings £j < kN/4, which does not intersect the perturbative strong 
tunneling regime. 

In the next section, we will discuss the influence of the gap on the quantum dynamics 
for the relative number of particles. In particular, we will show that the choice of couplings 
corresponding to the minimal and maximal points of the energy gap represent transition 
points separating different types of dynamical behaviour. Fortunately, we can identify all 
these points using eq. (jlj) above. 



3 Quantum dynamics 

We will investigate the quantum dynamics using the same methods we employed in [19]. 
The time evolution of any state is determined by {"^{t)) = U{t)\(f)o), where U is the tem- 
poral evolution operator given by U{t) = J2m=o |"^)("^| exp(— zi^^^), |"^) is an eigenstate 
with energy and |0o) represents the initial state. Using these expressions we can 
compute the expectation value of the relative number of particles 

((iVi - iV2)(t)) = {^{t)\N, - N^l^it)). (5) 
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In Fig. |3]we plot the expectation value for the relative number of particles in the Josephson 
regime {kjEj = 1) for different ratios of the coupling Afi/Sj for N = 50, 100 and with 
the initial state |A^, 0). 
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Figure 4: Time evolution of the expectation value of the relative number of particles in 
the Josephson regime k/Sj = 1 {Sj = 1), for N = 50, 100. The initial state is |A^, 0). The 
central graph corresponds to the configuration where the energy gap is a minimum, at 
Anc/Sj = 24.5, 49.5. The other graphs represent points before and after this minimum 
point given by Afi/£j = A^d^j T 5. 



The central graph in each column represents the dynamics when the system is in 
the configuration where the energy gap is minimal, which can be calculated through eq. 
(jD). This corresponds to Afi^/Sj = 24.5, 49.5 for the cases = 50, 100, respectively. 
The other graphs in each column, above and below, depict the dynamics for cases where 
Afi/Sj = Afic/Sj =F 5 in Fig. The coupling corresponding to the minimal energy gap 
signifies a point where the dynamics changes drastically: before it the behaviour is clearly 
non-periodic (the first box in each column), while after it the oscillations follow a collapse 
and revival sequence (the last box in each column). A noticible feature here is that in the 
same region, both graphs = 50, 100 show basically the same qualitative behaviour. The 
collapse and revival pattern, however, is more pronounced when the number of particles is 
larger. Further increasing the ratio Afi/Sj beyond the value corresponding to the minimal 
energy gap, the behaviour is still of a collapse and revival type, however the amplitude of 
oscillation and the collapse and revival time both decrease. 

Next we turn our attention to study the connection between the dynamics and the 
energy gap for the crossover region between the Josephson and Fock regimes, where a set 
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of minima and maxima in the energy gap appear as a function of the external potentiaL 
Again, the position of these minima and maxima can be obtained from eq. In Fig. 
El we plot the temporal evolution of the expectation value for the relative number of 
particles for the case = 50 for the coupling ratio k/Sj = 60, using the initial state 
|A^, 0). The central graph in each column represents the point where the energy gap is 
minimal (/ = 1,3) or maximal (/ = 2,4), occurring at A^JSj = 1470, 1440, 1410, 1380 
(from left to right). The other graphs above (respectively below) in each column represent 
the dynamics for the situation where Afi is smaller (respectively larger) than Afic- The 
first and third columns represent the dynamics for cases where the energy gap is near 
the two lowest minima, while the second and fourth represent the behavior near the two 
lowest maxima. 
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Figure 5: Dynamics of the expectation value of the relative number of particles for = 50, 
k/Sj = 60 {Sj = 1) and with the initial state |A^, 0). The central graph, at each column, 
represents the point where the energy gap is minimal (/ = 1,3) or maximal (/ = 2,4), 
occurring at Afic/Sj = 1470, 1440, 1410, 1380 (from left to right). In each case the graphs 
above and below show dynamics of points around the minimum or maximum, chosen as 
Afi/Sj = Afic/^j T 5. It is clear that each minimum and maximum signifies a change in 
the dynamical behaviour. 



The minima and maxima for the energy gaps represent frontiers between different 
types of oscillations in the neighbourhood of these points. After the last minimum at 
/ = 1 the expectation value evolves periodically, as in the last graph of the first column. 
However, the amplitude of the oscillations becomes smaller as the ratio Afx/Sj becomes 
larger. 
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Finally we investigate the dynamics for the situation where the system is below the 
threshold coupling {k/ £j < 4:/N) where the energy gap does not exhibit a local minimum, 
as already mentioned. In Fig. IHlwe plot the expectation value for the relative number of 
particles for different choices of the ratio Sj/k in the regime between Josephson and Rabi 
for = 40 and using the initial state |iV, 0). Here the dynamics displays collapse and 
revival of oscillations for any value of the external potential. The presence of an external 
potential increases the collapse and revival time, which tends to infinity as the external 
potential approaches infinity. Moreover, increasing the external potential also leads the 
system out of a delocalised phase into one of self-trapping. In fact, it is apparent that for 
any choice of Sj/k one can choose a suitably large A/i which will dominate the dynamics 
and the system will necessarily display self-trapping. 



ej/k=20 I N ^ 40 I ej/k=20 




Figure 6: Expectation value for the relative number of particles versus time for different 
values of the external potential in the regime between Josephson and Rabi (below the 
threshold coupling). Here,, we are fixing Sj = 20 and k = 1 The initial state is |A^, 0) 
with = 40. Increasing the external potential produces increasing collapse and revival 
time and also leads to a self-trapping phase. 



4 Classical dynamics 

Here we compare the above results with those that are obtained by an analogous classical 
analysis of the model. Defining z = {Ni — N2)/N we may equivalently consider the 
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Hamiltonian [26,30] 

H{z, 0) = ^ Q^' -Pz- cos(20/iV) j 



(6) 



where 



kN A/i 



and z, (f) are canonically conjugate variables. The Hamiltonian (j6j) obeys the symmetries 

In the usual way the dynamics is given by Hamilton's equations 

= ^ = ^(xz-(3+^=L=cos{2<p/N)] (7) 



dz 2 V VT 
OH 



z 



d(/) 



-Sj (^Vl - sin(20/A^) 



The classical dynamical evolution of the system is constrained to the level curves of H 
in phase space. The nature of the level curves is largely determined by the fixed points of 
H determined by i = = 0, and transitions in the dynamical behaviour can be identified 
with fixed point bifurcations. We find the following classification of the fixed points: 

• = and z is a solution of 

A^ = -^i= + /3 (9) 

which has a unique real solution for A > 0. For this solution the Hamiltonian attains 
a local minimum. 



(p = Nit/2 and z is a solution of 



This equation has either one or three real solutions. 



a) < A < 1 

For any value of f3 there is just one real solution, for which the Hamiltonian attains 
a local maximum. 



b) A > 1 

Here a transition value for Pc appears, which is dependent on A. For /? < Pc, the 
equation has two locally maximal solutions and one saddle point, while for (3 > j3c 
the equation has just one real solution, a locally maximal fixed point. Thus a fixed 
point bifurcation occurs at Pe- 
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Figure 7: Classical time evolution of the relative number of particles in the Josephson 
regime k/£j = 1, for = 50,100. The initial condition is 2; ~ 1, = 0. The central 
graph corresponds to the configuration where for the quantum model the energy gap is 
a minimum, at A/ic/^j = 24.5, 49.5. The other graphs represent points before and after 
this minimum point given by Afi/Sj = Afic/^j T 5. 



We remark that in the case A = 1 the transition value is given by (3c = 0. Thus in the 
absence of the external potential the classical dynamics predicts a transition coupling of 
k/Sj = 2/N [26], whereas for the quantum dynamics the threshold between delocalisation 
and self-trapping occurs at k/Sj = A/N [16,19]. The difference between the classical and 
quantum cases may be explained in terms of the uncertainty relation between z and (p. 

It is apparent that the analysis of the classical system for non-zero external potential 
only predicts one transition point, which can be identified as a fixed point bifurcation. In 
constrast we have shown above that for the quantum dynamics there are several points 
associated with minima and maxima of the energy gap, for which there is a qualitative 
transition in the dynamical behaviour. In Figs. (j7l8l9|) below, we show the dynamical 
behaviour predicted by the classical Hamiltonian using the same initial conditions and 
coupling parameters for the cases shown in Figs. (|4I5I6|) respectively. In particular, in 
figure ini we present both, the classical and the quantum curves in a shorter time interval 
for a better comparison. 

In most cases the results predicted by the quantum and classical analyses agree for 
sufficiently small time intervals, a feature also observed in [16]. In all cases however, there 
are stark differences in the dynamics over long time intervals. The irregular behaviour 
appearing below the transition coupling in Fig. |3]does not occur in Fig. [7| The dynamical 
transitions that are shown in the last three columns of Fig. El are not apparent from the 
classical results of Fig. |H1 And the collapse and revival of oscillations which is clear in the 
quantum dynamics (FigO and the solid curve in FigE)) is not seen in the classical case 
given by the dashed curve in Fig. |H1 
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Figure 8: Classical dynamics of the relative number of particles for N = 50, k/Sj = 60 
and with the initial condition 2; ~ 1, = 0. The central graph, at each column, represents 
the point where for the quantum model the energy gap is minimal (/ = 1, 3) or maximal 
{I = 2,4), occurring at An^/Sj = 1470, 1440, 1410, 1380 (from left to right). In each case 
the graphs above and below show dynamics of points around the minimum or maximum, 
chosen as A/i/Sj = An^/Sj =F 5. 



5 Summary 

To summarise, we have studied the nature of the energy gap as a function of the exter- 
nal potential and its connection to the dynamical behaviour. We have shown that the 
energy gap function presents no minimum below the threshold coupling, while above this 
threshold point one minimum occurs close to the Josephson regime and a set of minima 
and maxima appear in the Fock regime. We found that the total number of these minima 
and maxima depends only on the number of particles. An explicit expression for the 
position of these minima and maxima was given. The connection between these energy 
gap extrema and the dynamics for the relative number of particles was investigated, and 
it was found that the extrema determine transition points separating different dynamical 
behaviours. We have also compared these results with those that are obtained from a 
classical analysis of the Hamiltonian. 

Finally, we refer to recent experimental work [31] where both delocalisation, with 
tunneling times of the order 50 ms, and self-trapping have been confirmed in a single 
bosonic Josephson junction with a symmetric double-well potential using a Bose-Einstein 
condensate of ®^Rb atoms. In this experimental setup both scenarios were observed by 
simply varying the total atomic population A^. We note that from equation (jH) it can 
be seen that for fixed coupling parameters fc, A/i and Sj, the transition couplings A/ic 
can be tuned by also varying the total population A^. Thus the results we have obtained 
above are relevant for the experimental system of [31] with an asymmetric potential. 
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Figure 9: Classical (dashed line) and quantum (solid line) evolution of the relative number 
of particles versus time for different values of the external potential (from the top to the 
bottom: A/i/fc = 0,1,3 on the left and A/i/fc = 5,10,15 on the right) in the regime 
between Josephson and Rabi (below the threshold coupling). Here we are fixing = 
4:0,8 J = 20 and k = 1. The initial condition is z ~ 1, = 0. 
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